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INTRODUCTION 
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O T he J-point correlation functions have proved very useful 
I for the analysis of the large scale structure. One might 

P also use other approaches, such as the one-point probability 
• distribution function (PDF) (e.g. Gaztanaga & Yokoyama 

^ 1993, Sutherland, Maddox & Efstathiou 1994) but corre- 
lation functions have some clear advantages. For catalogs 
of angular positions the correlation analysis can be used to 
sort out the effects of projected clustering and estimate the 
intrinsic three-dimensional properties. This seems more dif- 
ficult (if not impossible) with the PDF. Moreover, there are 
analytical results for J-point correlations both in perturba- 
tion theory and in highly non-linear regime of gravitational 
growth (see below). Some of these models give degenerate 
solutions for J < 4, and they can only be tested using the 
higher order functions, e.g. equation (41) or (42). The J- 
point amplitudes can also be used as an expansion to the 
PDF (Juszkiewicz et al. 1993) with an accuracy given by 
the highest order J for which the amplitudes are known. 
Thus, there are important motivations to extend the anal- 
ysis of correlation functions to the largest orders available. 
The problem in practice is that for J > 3 the estimation of 
£j(ri, ...,rj) from galaxy catalogs is difficult as one has to 
work in a multidimensional space. In this paper we simplify 
the analysis by considering volume and area averaged corre- 
lations which gives better signal to noise properties and yet 
provides enough information to test the models. With this 



ABSTRACT 

We estimate J-point galaxy averaged correlation functions uJj(6) for J = 2,..., 9, in 
a sample of the APM Galaxy Survey with more than 1.3 x 10 6 galaxies and a depth 
T> ~ 400 Mpc. The hierarchical amplitudes sj = Uj/U^ -1 are roughly constant, 
up to J = 9, between 6T> ~ 0.5 /i -1 Mpc and 6T> ~ 2/i _1 Mpc and decrease slowly 
for larger scales. At scales 6T> > 7/i -1 Mpc we find strong similarities between the 
statistical properties of the galaxy fluctuations 6 g and the theoretical properties of 
matter fluctuations 6 m evolving under the influence of gravity in an expanding universe 
on assumption that the initial fluctuations are small and Gaussian. This is most easily 
explained if at large scales there is no significant biasing between matter and galaxy 
fluctuations, i.e. 6 g ~ 6 m . 

The comparison of the skewness in the Cf A and SSRS catalogues with comparable sub- 
samples of the APM indicates that the volume of a "fair sample" has to be much larger 
that the one in the combined CfA/SSRS catalogues. 
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simplification we are able to estimate up to the 9th order 
amplitudes of { j at small scales and up to the 6th order at 
larger scales. 

Observations of clustering in galaxy catalogues seem to 
find a particular relation or hierarchy for the J-point galaxy 
correlations, the so called Hierarchical model, which could 
be a trace of the gravitational nature of the evolution of 
large scale matter fluctuations. Early observational evidence 
for this hierarchy was found in angular catalogues of op- 
tically selected galaxies (e.g., Groth & Peebles 1977; Fry 
& Peebles 1978; Sharp et al. 1984). Recent analyses have 
confirmed these results and extended them to redshift cata- 
logues and to IRAS galaxies (e.g Szapudi, Szalay & Boschan 
1992; Meiksin, Szapudi & Szalay 1992; Gaztanaga 1992; 
Bouchet et al. 1993). From theory, hierarchical forms have 
been obtained in non-linear perturbation theory with ini- 
tial Gaussian fluctuations (Peebles 1980, Fry 1984b, Goroff 
et al. 1986, Bernardeau 1992) and for scale-free models in 
the highly non-linear regime of gravitational clustering (e.g. 
Davis & Peebles 1977, Fry 1984a, Hamilton 1988). Numeri- 
cal simulations also show similar properties (e.g. Efstathiou 
et al. 1988, Bouchet & Hernquist 1992, Weinberg & Cole 
1992, Lahav et al. 1993, Fry, Melott & Shandarin 1993, 
Lucchin et al. 1993). On the other hand, there is a large 
literature on models with non-Gaussian initial fluctuations 
which do not have hierarchical correlations at large scales 
(e.g., Moscardini, et al. 1993 and references therein). Such 
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models can display interesting behavior of the higher or- 
der moments, and observations might be able to distinguish 
them from models with initially Gaussian fluctuations. 

In the standard high peaks biasing model the observable 
galaxy distribution is related to the underlying matter dis- 
tribution assuming that galaxies form at peaks above some 
global threshold in the smoothed linear density field (e.g. 
Bardeen et al. 1986). In the limit of high threshold and small 
variance, this model is well approximated by the linear bias 
scheme, in which the smoothed galaxy 8 g (x) and mass den- 
sity fields 8 m (x) are linearly related: 8 g (x) = b 8 m (x). How- 
ever, since galaxy formation is a complex, non-linear process 
involving both gravitational and non-gravitational interac- 
tions, the relation between the mass and the galaxy distri- 
butions may be more complicated than in the high peaks 
biasing model. Fry & Gaztanaga (1993) have shown that 
any local biasing 8 g = f(8 m ) between the galaxy and the 
matter fluctuations preserves the hierarchical structure at 
large scales but changes the amplitudes. On the other hand, 
Frieman & Gaztanaga (1994) have shown how one could use 
the observed 3-point amplitude to identify the effects of non- 
local biasing. Thus the hierarchical amplitudes in the galaxy 
distribution can help us to learn about the corresponding 
matter amplitudes and the effects of biasing. 

In this paper we present a preliminary analysis of the 
the J-point correlations in the APM catalogue. Section §2 
contains the assumptions and the methods used in this pa- 
per while in section §3 we describe the APM sample and 
the results of the analysis. In section §4 we make a similar 
analysis over the CfA and SSRS catalogues to compare with 
the APM. The last section is devoted to a discussion of the 
implications. 



2 METHOD OF ESTIMATION 
2.1 Definitions and notation 

A simple way to define the J-point correlation function of 
a density distribution is by considering the density fluctu- 
ations. We define the density fluctuation 8 t = 8(r t ) at a 
point n by p t = p (1 + 8 t ), where p t = p(r t ) is the density 
at a point n and ~p is the averaged density: ~p =< p t >. 
The statistical average, denoted by < ... >, is over differ- 
ent realizations of 8(r) and corresponds to the average over 
positions in a fair sample of the universe. 

Using this notation the J-point correlation functions are 
defined as: 



bj >c 



(1) 



where < ... > c stands for the "connected" part of the ex- 
pectation values. The "connected" part corresponds to the 
contribution to the probability < 8\,...,8j > which does 
not include any conditional probability of lower order. Up 
to J = 3 we have < ... > = < ... > c but for J = 4: 



< 8 1 8 2 8 3 8 i >=< 8 1 8 2 8 3 8 i > c + > < 8,8j > < 8 k 8i > (2) 



E 

ijkl 



as the conditional probability to have a pair i, j given a pair 
k, I contributes directly to the probability < 81828384 >. 

For a point process the density n is defined by the prob- 



ability, dPi , that a galaxy is found to be inside a randomly 
placed volume element dV: 



dPt = n dV. 



(3) 



To relate the definitions in a continuous density field to the 
point process it is useful to consider the so called Poisson 
model, by which the probability dP(r) to find a galaxy in 
a given volume element dV around r is given by the local 
density, p(r), of a continuous field: 



dP(r) = p(r)dV, 



(4) 



and is independent of what happened in neighboring ele- 
ments. The probability to find one galaxy in a randomly 
placed element is dP\ =< dP(r) > and therefore, compar- 
ing equations (3) and (4), we have that n = ~p. This model 
gives a prescription to estimate the galaxy correlation func- 
tions. For example, the two-point correlation function, 6( r )' 
is given by the joint probability dP2 that two galaxies are 
found in the two volume elements dV\ and dV2 placed at 
separation r = T\2'- 

dP 2 = <pi P 2> dV 1 dV 2 = p 2 [ 1+ < 8^2 > ] dV 1 dV 2 

= n 2 [ 1 + 6(0 ] dV 1 dV 2 . (5) 

In a similar way, the three-point correlation function, 
6(^1 , f'2, f 3), is given by the joint probability dPs that three 
galaxies are found in the volume elements dV\ , dV2 and dVs : 



dP 3 = rf [ 1 + 6(12) +6(23) +6(13) 
+ 6(123) ] dV 1 dV 2 dV 3 , 



(6) 



where 6(u) = and 6(123) = 6(^1, r 2 , r 3 ). 

In the angular distribution the 2-point angular correla- 
tion function W2(0) is given in terms of the probability of 
finding two galaxies centered in each of the elements of solid 
angle cifii and d&2 at angular separation: 9 = 812: 



dV 2 = A/" 2 [ 1 + w 2 (0) ] rffii dtt 2 



(7) 



where M is the mean number of galaxies per unit solid angle. 
For the 3-point angular correlation, (113(^2,^3), 

dV 3 = AT 3 [ l + w 2 (12) + w 2 (23) + w 2 (13) 

+ w 3 (123) ] d^d^dtts- (8) 

Using correlation functions is particularly useful for angu- 
lar distributions as they conveniently sort out projection 
effects. This can be illustrated for example in the above 
equation. The first term in equation (8) accounts for uncor- 
rected triplets that are clustered together just by chance or 
because of the projection. The next three terms correspond 
to a combination of a correlated pair with an uncorrelated 
galaxy that forms a triplet also by chance or projections, 
whereas the last term corresponds to real clustering. Thus, 
by studying 103 instead of dV 3 we avoid projection effects. 
Consider now the (smoothed) density contrast 8w(x): 



8w(x) 



1 

Vw 



dx'S(x')W(x 



dx W(x) 



(9) 



(10) 



T/U(a;) is the window function and the integral is over all 
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space <S. For a top-hat window, 8w(x) is the volume average 
of 8(x). In terms of the smoothed fluctuations, the volume- 
averaged J-point correlation functions are: 



tj(Vw) = < hw >c 



(11) 



Because of the fair sample hypothesis < ... > corresponds to 
the average over positions in the sample. Therefore, { j does 
not depend on x and it is just a function of the window and, 
in particular, of the volume Vw in equation (10). It follows 
from equation (9) that: 

J J = -L I dr 1 ...drjW(r 1 )...W(rj)ij(ri,-,rj). (12) 
v w J s 

The case J = 2 is commonly used to characterize rms fluc- 
tuations = (, 2 (Vw), 



«2(^v) = ^y dkkP(k)w 2 (k) (is) 

in terms of the power spectrum of fluctuations P(k). 



2.2 Area-averaged correlations 

The difficulty of estimating J-point angular correlations, 
wj(9i,...,9j), for J > 2, and the resulting amount of in- 
formation are overwhelming as one has to consider all pos- 
sible configurations in a multidimensional space. We do not 
necessarily need all this information. A considerable simplifi- 
cation can be achieved by estimating the hierarchical ampli- 
tudes from area-averaged correlations Toj, which simplifies 
the data analysis and gives better signal to noise properties. 

We define the area-averaged angular correlations TJj(9) 
in terms of the angular correlation functions wj(9i, ...,9j): 



i 7 



dAu.JAj wj(9 1 ,...,9j) 



< 8 J (6) > c 



(14) 



where A = 27r(f — cos#) is the solid angle of the cone, dAj = 
sin jdO jdip j and 8(9) are the fluctuations inside the cone. 
Thus To ' j(9) only depend on the size of the cone, 9, as they 
correspond to smoothed correlations. In terms of the spatial 
correlation functions, £,j(ri , r j), 

Toj(9) = / dr 1 ...drjip(r 1 )...ip(rj)£j(r 1 ,...,rj), (15) 
Jv 

where dr t is the proper volume element (see §2.4), ij)(r) i s 
the normalized probability that a galaxy at a distance r is 
included in the catalogue, i.e. equation (27), and V is a cone 
of angle 9 and infinite depth. 

In our analysis we estimate Toj(9) from the central mo- 
ments of the angular counts 



TO J 



(9) = J^(i-N) J Pi(9), 



(16) 



where Pi (6) is the probability of finding i galaxies in a ran- 
domly selected cell of solid angle A = 27r(f — cos#) and 
N = i Pt. Note that the galaxy density fluctuation in 
the cell is 8 g = (i — N)/N, and therefore 



TOJ = N j <8g > 



To obtain the correlation functions To j we have to estimate 
the "connected" moments fij = N < 8g > c . The relations 
between raj and (ij up to J = 9 are given in the Appendix 
Af . Because of the discreteness < 8 J g > c is not a good esti- 
mator of To j unless N ^ I . A better estimator of To j is given 
in the Appendix Af. 



2.3 Hierarchical amplitudes 



In the hierarchical model all high-order correlations can be 
expressed in terms of the two-point correlation function: 



,rj) 



j-i 

ab 



(18) 



In graphical notation, associated with each term in equa- 
tion (18) there is a graph, such that vertices, or nodes, cor- 
respond to the points T\ , . . . , r j, and edges, or lines, between 
node a and node b correspond to factors {2 (fab) that connect 
all points. Thus the hierarchy is composed of "tree" graphs 
(connected with no cycles) of J vertices and ,1 — 1 edges. 
The sum over a denotes topologically distinct graphs; the 
sum over ab is over relabelings within a. At each order J 
there are in total J J ~ 2 terms, corresponding to all possible 
reassignments of the labels a, b = 1, . . . , J. Sometimes the 
labels a in Qj, a are omitted and Qj is the average over 
different topologies: 



Qj 



J 



(19) 



where N a is the number of graphs with topology a. 

For the average correlations fj(-R) in a spherical cell of 
volume V = 4/3wR 3 , i.e. equation (12), 



fj(-R) = — / dr 1 ...dr J ij{r\, ...,r j) 



V J 

the hierarchy (18) translates into: 

IAR) = Sj [I 2 (R)] j -\ 

where Sj are related to Qj by: 
Sj = BjJ j - 2 Qj. 



(20) 



(21) 



(22) 



Parameters Bj are estimated in Appendix A2. For the area- 
averaged correlations (15) we define: 

TJj(9) 



sj(9) 



[TJ 2 (9)y-i- 



(23) 



(17) 



The relations between sj and Qj or Sj follow after a careful 
study of projection effects. 



2.4 Projection effects & selection function 



To estimate the spatial properties of density fluctuations we 
use angular positions of galaxies over a "complete" magni- 
tude limited survey. The analysis takes into account two dif- 
ferent projection effects. One is to distinguish between real 
clustering and clustering from projections. Angular correla- 
tion functions conveniently sort out this problem (see section 
§2.1). The other projection effect comes from the fact that 
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Table 1. Projection factors for different slopes 7 and parameters in the luminosity function: M£ and ag. 
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catalogues are limited by apparent magnitude which intro- 
duces a selection function. 

The metric for an homogeneous universe, in comoving 
coordinates x, is: 



dr 2 = c 2 dt 2 — a(t) 2 [f 2 (x)dx 2 + x 2 dA 2 ] 



(24) 



where dA 2 = dd 2 + sin 2 9 dip 2 . The expansion factor a = 
(1 + is given in terms of the redshift z = z(x) which 

follows Hubble's law: 

H x/c = \q z - go + 1 + (go - + 2q z\ / (1 + z), (25) 



where go = fio /2 in a matter-dominated universe. The Hub- 
ble parameter Ho and all units are given in terms of h so 
that Ho = WOh km s" 1 Mpc" 1 . The proper volume element 
is dr = a 3 _F _1 (x) 



F(x) = [1 - (H x/cf(Q - 1)] 



1/2 



(26) 



is the correction for curvature. 

The selection function ip( x ) i s the normalized probabil- 
ity that a galaxy at coordinate x is included in the cata- 
logue. This probability depends on x because given a galaxy 
with absolute magnitude M its apparent magnitude m is 
a function of x. We make the usual assumption that there 
exists a universal luminosity function (f>(q), so that absolute 
magnitudes and positions of galaxies are uncorrelated. For 
a catalogue with apparent magnitudes between mi and rri2 
the selection function is: 



ip(x) 

q t (x) 
Mi{x) 



dq 4>{q) 



10 -f(M, W -M*) i=lj2 

ib, — 51og 10 di(x) — 25 

/ * a —q 

<p q e . 



(27) 



where di = x(l + z) is the luminosity distance in ft -1 Mpc 
and 4>(q) is the Schechter luminosity function. We choose 
the amplitude <j>* to be a constant (the value of which turns 
out to be irrelevant in our analysis) whereas M* = M*(z) 
and a = a(z) are a function of redshift z to account for k- 
corrections and evolution. Following Maddox et al. (1990a) 
we use: 



M* = Mo + M* 



2.5 Projections and hierarchical amplitudes 

Because of the projection and the selection function ip( x ) 
there is a complex relation between the angular amplitudes 
sj in equation (23) and the corresponding Sj or Qj in three 
dimensions. Nevertheless, given a model for the two-point 
correlation function {2, it is possible to find a numerical 
estimate for sj. The contribution from the angular average 
(14) is very small whereas the one from the projection can 
be estimated as in Groth & Peebles (1978). For the hierarchy 
(18) with a power-law correlation ^(t) = (ro/r) 7 and small 
angles, 9 <C 1 in radians, 



sj 



rj Cj J J ~ 2 Qj 



(29) 



Parameters Cj ~ 1 correspond to the angular average in 
equation (14) and are estimated in Appendix A2, whereas 
rj are related to the selection function ip: 



rj 



h 



(30) 



F x 2 dx f k xt 3 -^- 1 ) (1 + 



a + oi\ z 



(28) 



The proper clustering ^(t, t) has been parametrized as 
F^2{t)(1 + 2) -( - 3+e -' where F = F(x) is the correction for 
curvature in equation (26). We use the selection function 
ij) = ip(x) given by equation (27) with the parameters by 
Maddox et al. (1990a), i.e. e = and M * = -19.8, M* = 1, 
ceo = —1, «i = —2 in equation (28). The resulting values 
of t j increase with 7 and M * and decrease with «o, but do 
not change much within the uncertainties in the shape of 
the luminosity function (see also Peebles 1980, §56). This 
is illustrated in Table 1 where values of rj are plotted for 
different parameters in the selection function. For a fixed 
7 = 1.7, the values of rj only change by a few percent even 
when Mo* and cto are changed by AM " = 1.0 and A«o = 0.4 
(the values in Table 1 correspond to the combinations that 
give the largest changes in rj). The ranges AM " = 1.0 and 
A«o = 0.4 are much larger than the observational uncer- 
tainties in the luminosity function (e.g. Efstathiou, Ellis & 
Peterson 1988) and have been exaggerated to illustrate how 
rj are insensitive to the detailed shape of 0. Small variations 
in fio = 1, e = 0, M* = 1 and a\ = —2 produce even smaller 
changes. It is important to notice that although rj are un- 
affected by changes in ip, the overall normalization oi Ik can 
change quite a lot. While the amplitude of {2 is uncertain 
by 40% for AM * = 1.0 and Aa = 0.4 the corresponding 
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uncertainty in T3 is only 2%. This is an excellent motivation 
for using the values of Sj as measures of clustering. 

In the analysis of the APM below, variations of 7 are 
only important for very large scales, 9 > 3°, where 7 changes 
from 1.8 to 3.. In this case r j suffers a considerable variation, 
see Table 1, and equation (29) is not a good approximation. 



3 ANALYSIS OF THE APM DATA 
3.1 The APM Galaxy Survey 

Galaxies from the APM Galaxy Survey (Maddox et al. 
1990a-c) are selected with apparent magnitudes, ra = bj, be- 
tween mi = 17 and m,2 = 20. We use equal area projection 
pixel maps with a resolution of 3.52 arcminutes containing 
over 1.3 x 10 6 galaxies. With the selection function above 
the depth, defined as the distance where the distribution of 
galaxy distances peaks, is D ~ 380 ft -1 Mpc. 

Residual contamination from merged images is / ~ 5% 
for our magnitude range (Maddox et al. 1990a). Assuming 
that these mergers are uncorrelated, their net effect is to 
dilute the clustering. The observed fluctuations, 8 , are: 



So = (1 " /) Sg + fS„ 



(31) 



where 8 g are the galaxy fluctuations and 8 m the fluctuations 
from the mergers. Because the mergers are uncorrelated < 



5^ > = < 8^8f >= for all N, M, and therefore: 



<8 J > C = (l-f) J <6l>c 



(32) 



which means that the observed values of the hierarchical 
amplitudes s°j =< 8„ >c / < 8 2 >c -1 have to be multiplied 
by (1 — /) ~ 2 to obtain the galaxy amplitudes, i.e. 



sj= (I-/) 7 " 2 S°j 



(33) 



These factors are approximate and therefore the error-bars 
will be at least as large as the corrections. 

3.2 Counts-in-cells and error estimation 

Counts-in-cells are estimated for square cells of side I in a 
range I = 0.03° — 20°. The smaller scale is chosen to avoid 
shot-noise fluctuations. Starting from a cell where N ~ 1 
we select different cell sizes so that the area of the next cell 
is at least twice as big as the previous one. To estimate the 
errors we try two prescriptions. 

The first prescription is to estimate the dispersion in 4 
different zones within the sample. In this case, the number of 
test-cells we have used to obtain the counts P t is arbitrarily 
large. This approach is very conservative as the errors cor- 
respond to the cosmic variance in 4 smaller surveys which 
will be significantly larger than the cosmic variance for the 
full survey. 

For the second prescription, we first do a partition of 
the map into cells of a given size. We then group the cells 
(at random) in subsets and estimate the variance from the 
results in each of the subsets. This prescription provides a 
way to estimate the sampling errors without introducing the 
cosmic variance of the smaller maps. We have checked that 



the number of subsets used to estimate the errors is not very 
important; similar results are found using 4, 8 or 20 subsets. 
This is because we are in fact estimating numerically the 
internal variance in each quantity. This estimation can also 
be done using the central moments raj. For example, the 
second moment &2 = To~2N in equation (A6) can be written 
as: 



N) — N 



(34) 



where the sum extend over the M cells and N t is the number 
of galaxies in cell i. If the cells are independent, the variance 
of &2 is: 



Var(k 2 ) =< k\ > - < k 2 > 2 



2wi3 + m.2 — ra\ 



M 



(35) 



where raj are the central moments in equation (16). For 
higher order moments, kj, the corresponding expressions 
become more complicate and depend on higher orders raji, 
J' = 2, 2 J. In practice it is better to estimate numerically 
the variance by using the variance of the subsets. We have 
checked that for the lower moments J < 5 both estimations 
agree well. 

About 3% of the pixels correspond to "small holes" 
caused by unmatched images (satellite trails, bright stars 
or very bright galaxies) which are more or less randomly 
distributed. To estimate fluctuations we accept a test-cell 
provided the holes inside cover less than ~ 30% of the cell 
and compensate the deficit by 'filling' the hole with the lo- 
cal density. This assumes that the clustering of the missing 
part of the cell is identical to clustering in the rest of the cell, 
which is a rough approximation as different scales may have 
different clustering. Nevertheless, this approximation does 
not affect the results within the resolution of our analysis 
because each cell is a 100% bigger than the previous one. 



3.3 Two-point correlations 

We first estimate the two-point angular correlation function 
W2(0) in equation (7) from counts in the pixel maps using: 



w 2 (6) 



(36) 



that is, given a cell i with N t galaxies we find all cells j sep- 
arated by an angle 6 + dd and estimate the number of corre- 
lated pairs. The result is shown in Figure 1, which also shows 
the average correlation functions LO2(0) estimated from mo- 
ments of counts-in-cells. Error-bars correspond to the dis- 
persion in T02 estimated from narrow stripes that divide the 
sample in 4 zones (same zones as used by Baugh & Efs- 
tathiou 1994). The values of 102(d) are plotted as a function 
of the radius 6 = lj\pK of a circle with the area of cell. Fig- 
ure 1 also shows as a continuous line lo%(8) estimated from 
a numerical integration of W2(0), i.e. using equation (14). 
The agreement between the two estimations indicates that 
square cells give very similar results to circular cells once 
the size of the cells are scaled to 6 = l/i/w. For clarity, we 
have not shown the errors in the estimation of W2(0) which 
are quite big at large scales (see Maddox et al. 1990a). The 
"break" from the power-law is more gentle in T02 than in 
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lO- 3 - 
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Figure 1. The two-point correlation ~2(0), squares with error- 
bars, estimated from counts-in-cells in the APM map compared 
with W2(0), dashed line. The continuous line correspond to ~|(^) 
estimated from a numerical integration of W2 (0). 



%»2 because the averaging in equation (14) mixes small with 
large scales. 

A fit of the two-point angular correlation to a power-law 



W 2 



AO 1 ' 1 , for scales 9 < 2° gives A~2.7x 10" 2 and 7 



1.7. Using the inversion presented in §2.5 the corresponding 
two-point correlation function is {2 — (/o /r) 1 with ro ~ 
5 ft -1 Mpc. As mentioned above the largest uncertainty in 
this value is the one from the selection function. 



3.4 Three-point averaged correlations 

The three-point averaged correlation 103(8) is plotted in Fig- 
ure 2. Also shown is the hierarchical amplitude S3 = Us/Zu^j 
which is divided by 10 2 to be on the same scale. The er- 
rors are obtained from the variance in 4 different zones. The 
dashed line is the mean of the values in each zone, while the 
triangles correspond to the full map. At scales 9 > 2°, the 
mean of the zones deviates from the value in the full map 
which indicates that the cosmic variance in the smaller zones 
is important at these scales. As a consequence the errors at 
larger scales are too conservative. 

In what follows we use our second prescription for the 
errors that estimates the variance in 4 subsets of the parti- 
tion of the map. This prescription gives smaller errorbars at 
large scales. The new estimation is shown in Figure 4 below. 

Figure 2 shows that 33(8) is roughly constant for a large 
range of scales; 33(8) only changes a factor ~ 2 when we 
change 9 by a factor of ~ 10 2 . For the smallest scales in 
Figure 2, i.e. 9 < 0.1°, the amplitude S3 seems to turn over 
and decrease again but this effect is probably caused by 
shot-noise (see Appendix Al). 

The three dimensional amplitude S3(R), with R = 9V, 
corresponding to 33(9) is given by: 

33(9) (1-/) 5 3 (7) 



10" 1 



io- 2 



10- 




Figure 2. The three-point correlation "3 (6), squares, compared 
with the hierarchical amplitude S3 (6) = ZJ3 f~\ , triangles, divided 
by 10 2 to be on scale. 




m 3 



S 3 (9V) 



T3(j) C 3 (j) 



(37) 



R/8D (Mpc h- 1 ) 



Figure 3. Amplitudes 53 (R) for spherical cells of radius R = 8V, 
estimated from £3(0) in the APM. The dashed line corresponds 
to the amplitude S3 expected for the APM in the gravitational 
instability scenario using second-order perturbation theory. 



from equation (29) and equation (22), where (1 — /) is the 
merging correction to S3, i.e equation (33). This expression 
is accurate for 9 < 2° as the power-law model in {2 is a good 
approximation (see §2.5) and the hierarchical amplitude is 
roughly constant. For larger scales, 9 > 2°, we approximate 
53 with equation (37) by using a local slope 7 = j(9) (from 
Figure 2). The uncertainties in S3(#) at these large scales 
are larger or comparable to the error in this approximation. 
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3.4-1 Comparison with perturbation theory 

In order to compare these estimations with the values 
in non-linear perturbation theory, we calculate Ss(R) = 
£ 3 (i?)/£ 2 (-S) with f 2 given by equation (13) and: 



3 



d 3 k 1( i 3 k 2 B(k l ,k 2 ,\k l +k 2 | 



W(hR) W{k 2 R) W(|ki +k 2 |i?) 



(38) 



(see also Frieman & Gaztanaga 1994). For the bispectrum, 
5i23 = B(ki, k 2 , IC3), we use the second-order perturbation 
theory result: 



B12 



10 

T + 



ki ■ k-2 
k\k 2 
k - 2 



tl _|_ 

k 2 k\ 



k\k 2 



PiP 2 + (1 ~ 2) + (2 ~ 3) (39) 



(Fry 1984), where the power spectrum P t = P(k t ) is the 
one in linear perturbation theory. We estimate S3 numer- 
ically using the values of P(k) obtained by Baugh & Efs- 
tathiou (1993, 1994) for the same APM sample. Of course, 
this measured power spectrum corresponds to the non-linear 
evolution, but at large scales, R > To, it provides a good ap- 
proximation to the linear P(k). The resulting S3(R) is plot- 
ted as a dashed line in Figure 3 where it is compared with 
the APM skewness. There is a good agreement for scales 

5 h' 1 Mpc < 9V < 30 h' 1 Mpc. The last point in Figure 3, 
at 9T> ~ 40 ft -1 Mpc, is also consistent with the prediction 
after allowing for the large uncertainties in P(k) at large 
scales, which are not plotted (see Baugh & Efstathiou 1993 

6 1994). At small scales, 9V < 5 ft -1 Mpc, the galaxy skew- 
ness is significantly bigger than the perturbation result. 

The agreement in Figure 3 demonstrates that the clus- 
tering in the APM is consistent with the idea that fluctua- 
tions grow according to the gravitational instability picture. 
The only other assumptions that has been used in this com- 
parison, i.e. in equation (39), is that the initial conditions 
were Gaussian. 



3.5 Higher order correlations 

Correlations Zoj(9) up to J 



terms of the hierarchical amplitudes sj = toj/toi 1 . The er 



9 are shown in Figure 4 in 
j = Tc jjTc 3 2 ~ 

rors are from 4 subsets of the map. Higher order correlations 
J > 9 are dominated by poor statistics. 

The hierarchical amplitudes sj in Figure 4 are roughly 
constant at small scales 9 ~ 0.2° and decrease slowly for 
larger scale. For the smallest scale in Figure 4, 9 < 0.05°, 
the amplitudes seem to turn over and decrease again but 
this effect is probably caused by shot-noise (see Appendix). 

Table 2 shows the averaged values of the amplitudes be- 
tween 0.05° < 9 < 0.3°, i.e. 0.3 h' 1 Mpc < 9V < 2 h' 1 Mpc. 
The values for s J j FM are obtained from Figure 4 and the 
corresponding sf FM are compensated from merging and 
projection effects using: 

fV' 2 Bj(y) 



Sj(9V) 



sj(9) (1 



(40) 



from equations (29) and (22), with (1 — f) J ~ 2 the merg- 
ing correction to sj, i.e. equation (33). These relations are 



10 4 



10 E 



10 2 



10 1 



nrm 1 — 



I I. 



n 



1 




0.1 1 

8 (deg.) 



0.1 1 
9 (deg.) 



Figure 4. The hierarchical amplitudes sj = ~j/~ 2 ~ for the 
APM sample. Error-bars correspond to the dispersion in 4 differ- 
ent zones in the APM map. 



Table 2. The hierarchical amplitudes in the APM averaged be- 
tween 0.05° < 6 < 0.3°, i.e. 0.3 h' 1 Mpc < 6V < 2/i _1 Mpc. 



„APM 
J 



gAPM 



-)APM 



3 


4.95 ± 0.10 


4.04 ± 0.08 


1.31 ± 0.03 


4 


59.2 ± 2.4 


37.3 ± 1.5 


2.12 ± 0.09 


5 


1081 ± 59 


502 ± 28 


3.5 ± 0.2 


6 


(2.5 ± 0.1) X 10 4 


(8.7 ± 0.7) X 10 3 


5.4 ± 0.4 


7 


(6.6 ± 0.4 X 10 5 


(1.7 ± 0.2) X 10 5 


7.4 ±0.7 


8 


(1.8 ± 0.1) X 10 7 


(3.2 ± 0.4) X 10 6 


8.4 ± 1.0 


9 


(4.4 ± 0.4) X 10 8 


(5.8 ± 0.8) X 10 7 


7.7 ± 1.1 



only valid assuming power-law correlations and constant Q j . 
For these scales the power-law model in the two-point cor- 
relation is a very good approximation and the hierarchical 
amplitudes can be well approximated by constant values as 
shown in Figure 4. The errors in Sj and Qj are the max- 
imum between the errors from sj and the uncertainties in 
the merging correction. 



Table 3. Hierarchical amplitudes in the APM averaged at large 



> 1° 
J 


, i.e. V6 > 7h- 

„APM 
J 


- 1 Mpc. 

gAPM 


Qj PM 


3 


3.81 ± 0.07 


3.16 ± 0.14 


1.02 ± 0.05 


4 


32.5 ± 4.2 


20.6 ± 2.6 


1.17± 0.15 


5 


384 ± 62 


180 ± 34 


1.25 ± 0.24 


6 


3260 ± 1340 


1150 ± 580 


0.71 ± 0.37 
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if 



10 5 =- 



1 10 
R (Mpc h" 1 ) 



1 

R (Mpc h" 1 ) 



Figure 5. The hierarchical amplitudes Sj = inverted 
from the angular amplitudes in Figure 4. 



Table 3 shows the best fit of the amplitudes in Figure 4 
to a constant, weighted by the errors, for scales 9 > 1°, i.e. 
T>8 > 7 ft -1 Mpc. Notice that the number of points and the 
range of scales are different for each J. For J = 3 there are 
8 points extending from 7 ft -1 Mpc to 80 ft -1 Mpc whereas 
for J = 6 there are only 2 points, between 7 ft -1 Mpc to 
10 ft -1 Mpc. We have only use the values where the errors 
are not compatible with zero. The amplitudes Sj and Qj 
are compensated for merging and projection effects using 
equation (40) and equation (22) with the values of rj, Bj 
and Cj corresponding to the local slope 7. 

Figure 5 shows the estimated values of Sj from equa- 
tion (40) to illustrate Tables 2-3. By approximating Qj to 
be constant in all the range we have probably introduced a 
small distortions in the shape of Sj between 2 ft -1 Mpc < 
V9 < 7 A -1 Mpc, where there is a small variation of the 
amplitudes with scale. 



3.5.1 Comparison with the Lick map 

The results in Tables 2-3 can be compared with the averaged 
values of Q j from the Lick map (Shane and Wirtanen 1967). 

i 1.25 ± 0.18 while 



Lick 



Groth & Peebles (1977) found Q^' c 
Szapudi, Szalay & Boschan (1992) estimated Q, 
0.07. For J = 4, Fry & Peebles (1978) found Qf" 



1.44± 
t 3.0 ± 



0.5 while Szapudi, Szalay & Boschan (1992) quote Q 4 lc ~ 
1.95 ±0.33. These results correspond to the values averaged 
for all scales smaller than ~ 5 ft -1 Mpc. But, according to 
our analysis, the amplitudes decrease for V9 > 2 ft -1 Mpc 
and averaging could be misleading. This might explain the 
discrepancies between the two estimations of Q\ ,ck . For J > 
4 the amplitudes in the Lick map found by Szapudi, Szalay 
& Boschan (1992) are quite uncertain, but up to J = 6 they 
seem compatible with the APM amplitudes at larger scales. 



3.5.2 Comparison with perturbation theory 

We are going to compare the amplitudes in Table 3 with 
the predictions in non-linear perturbation theory, as the re- 
sults from N-body simulations (e.g. , Efstathiou et al. 1988, 
Bouchet & Hernquist 1992, Weinberg & Cole 1992, Lahav 
et al. 1993, Fry, Melott & Shandarin 1993, Lucchin et al. 
1993) prove that the perturbative approach works well even 
for scales where {2 is only slightly smaller than unity, i.e. for 
scales larger than ro ~ 5 ft -1 Mpc. 

In linear perturbation theory for gravitational growth 
of initially Gaussian fluctuations, the dynamics are com- 
pletely described by the two-point correlation {2 or its power 
spectrum, P(k), which grow with the scale factor. Higher 
order correlations are zero in linear theory as the distri- 
bution keeps Gaussian. In non-linear perturbation theory 
higher order correlations follow the hierarchical pattern, 
in equation (21), so that at each stage of the evolution 



Sj 



in a self similar way. The matter ampli- 



tudes Sj t m are characteristic of gravitational instability in 
an expanding universe and have been calculated analytically 
by Fry 1984b, Goroff et al. 1986 and Bernardeau 1992. In 
the quasi-linear regime Sj tTn do not vary much with time or 
fio ( Juszkiewicz, Bouchet & Colombi 1993) and therefore all 
dynamics are completely described by the power spectrum, 
just as in the linear perturbation case. Thus, Sj tTn can be 
used to check if fluctuations grow via gravity. 

The values of Sj tTn in non-linear perturbation theory 
do depend slightly on the shape (but not the amplitude) of 
P(k) and therefore on the scale. This has been illustrated 
in Figure 3 where 53 >m was estimated from the measured 
P(k). At a given scale Sj tTn can be approximated if we 
know the local slope n in P(k) ~ k n . For quasi-linear scales 
the observed local slope is n ~ —1, which corresponds to 
7 ~ 2. Thus, we can compare the amplitudes in the APM, 
Sf FM with the ones in perturbation theory 5j >m expected 
for P(k) ~ A; -1 . Juszkiewicz, Bouchet & Colombi (1993) 
found 53 >m = 34/7 — (n + 3) for a spherical top-hat window, 
so that S 3 ,m ~ 2.9 for n ~ -1, close to Sf PM ~ 3.16 ± 0.14 
in Table 3. A detailed comparison using the exact prediction 
is shown in Figure 3. Lucchin et al. (1993) have found, using 
numerical simulations, S4 m ~ 16 — 18 and S5 m ~ 135 — 147 



for n 



-1, which are also similar to 5 4 



20.6 ± 2.6 



and S£ PM ~ 180 ± 34 in Table 3. 



Thus the amplitudes in the APM, at least up to J = 5, 
agree quite well with the predictions in the gravitational 
instability scenario. 



4 COMPARISON WITH THE CFA AND SSRS 



Values S( 



4 and S( 



20 - 40 in tables 2- 



3 are larger than the ones in the CfA and SSRS redshift 
catalogues: S3 ~ 2 and S4 ~ 5 (Gaztanaga 1992). This is 
probably not due to redshift distortions (Fry & Gaztanaga 
1994). Is this then an intrinsic discrepancy or is it caused by 
the use of different techniques, e.g., the selection function? 

To answer this we estimate S3 for the North Zwicky 
Center for Astrophysics Survey (Huchra et al. 1983, here- 
after CfA) and the Southern Sky Redshift Survey (da Costa 
et al. 1988 hereafter SSRS) using the same techniques as 
for the APM. In the CfA we select 1930 galaxies in a solid 
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DISCUSSION 



m 2 



n 1 1 — i — r 




SSRS * 
CfA § 

J L_ 



1 



BD (Mpc h" 1 ) 



10 



Figure 6. The hierarchical amplitudes S3 for APM sample, con- 
tinuous line, compared with the ones in the CfA sample, open 
squares, and the SSRS sample, closed squares. The dashed lines 
enclose a 70% interval of confident in the dispersion of the results 
of S3 in 32 independent zones inside the APM, each extending 
over a scale of 6V ~ 100 h' 1 Mpc. 



angle of 1.83 sr. with m B < 14.5, 8 > and b > 40°. In the 
SSRS, 1760 galaxies are selected in a solid angle of 1.75 sr. 
with logD(O) > 0.1, 8 < -17.5° and b 11 < -30°. We use 
circular cells for the counts but this should not make much 
difference as explained above. 

Results are presented in Figure 6 in terms of Ss^dV) 
using equation (37) with T>cfA = 40 ft -1 Mpc and Vssrs = 
50 A" 1 Mpc, in agreement with the comparison by da Costa 
et al. (1988). According to the luminosity function in the 
CfA (Efstathiou, Ellis & Peterson 1988) the projection fac- 
tor T3 ~ 1.4. Relation (37) provides a very good approxima- 
tion as {2 follows the power-law model for the scales we are 
interested, i.e. 8T> < 10 ft -1 Mpc, and the three-point ampli- 
tudes seem roughly constant. Errors in the CfA and SSRS 
are from 4 realizations of the cell positions. Within the er- 
rors, the CfA and SSRS results seem to agree well with the 
results of Gaztanaga (1992) and they have lower amplitudes 
than the APM. 

A natural explanation for the discrepancies between the 
APM and the CfA/SSRS catalogue is sampling effects, i.e. 
that our local neighborhood is not a fair sample. After all, 
the APM map covers over 200 times more volume and 300 
times more galaxies than the combined CfA and SSRS. 

To investigate this point we divide the APM map into 
a 8 x 4 grid, i.e. 32 zones each extending over a scale of 
V9 ~ 100 h' 1 Mpc, about the combined depths of the CfA 
and SSRS. We have compared the values of S3 in each inde- 
pendent APM zone to estimate a 70% interval of confidence 
in the values of S3. This is displayed in Figure 6 and com- 
pared with the CfA and SSRS results, which lay within the 
lower end of the APM interval of confidence. 



The hierarchical amplitudes sj in Figure 4 are roughly con- 
stant at small scales 9V ~ 1 ft -1 Mpc. This result agrees 
with the scaling relations expected in the similarity solu- 
tions of the BBGKY model for the highly non-linear regime 
(e.g. Davis & Peebles 1977). The values of Qj have been 
estimated for the similarity solutions with different approx- 
imations. Fry (1984a) found 



Qj 



J 



2,1-2 



whereas Hamilton (1988) estimated: 



Qj 



J ~ 2 J! 
2 ' 



(41) 



(42) 



Thus, to check these models we have to assume a value of 
Q3 and compared higher order amplitudes, J > 3. Even 
after allowing for large uncertainties in the value of Q3, the 
predictions in equations (41) or (42) are much smaller than 
the values of Qj FM in the APM, i.e. in Table 2. One might 
argue that these discrepancies could be caused by small scale 
biasing (see below). An interesting feature shown in Table 2 
is that the values of Q J j FM at small scales increase with the 
order J up to J = 9. 

Sutherland, Maddox & Efstathiou (1994) have found 
that at large scales the angular counts in the APM Survey 
fit the lognormal distribution. The lognormal distribution 
reproduces the Kirkwood scaling (see e.g. Coles & Jones 
1991) and therefore has amplitudes Qj that increase with 
{2 > 1. At large scales, {2 < 1, the amplitudes in the log- 
normal distribution converge to constant values Qj ~ 1, 
which are similar to the APM amplitudes, Qj PM , in Table 
3. According to perturbation theory there is nothing special 
about Q j ~ 1, as different values of the power index n yield 
different values of Qj. For n ~ — 1 we have that Qj ~ 1 
is a good approximation, but at larger scales where n > 0, 
Qj are significantly lower. Thus the agreement between the 
lognormal distribution and perturbation theory only occurs 
for the small range of scales where n ~ — 1. In any case, the 
amplitudes grow much too rapidly in the lognormal distri- 
bution and at 9V ~ 1 ft -1 Mpc, Q j are much larger than 
values we found in the APM sample. 

In section §3.3 we have shown that the amplitudes 
q apm ^ n fp^ig 3^ a veraged at large scales, agree quite well 
with non-linear perturbation theory. This results agrees with 
the comparison of S3 with second-order perturbation the- 
ory in Figure 3. There are therefore strong similarities be- 
tween the properties of galaxy fluctuations and the prop- 
erties of matter fluctuations evolving under the influence 
of gravity in an expanding universe. This indicates that 
the galaxy distribution has similar properties to the mat- 
ter distribution: i.e. the effect of biasing on Sj >m is small. 
Fry & Gaztanaga (1993) have shown that any local bias- 
ing 8 g = f(8 m ) = ^2bk8 k between the galaxy fluctuations 
8 g and the matter fluctuations 8 m preserves the hierarchical 



structure at large scales but changes the amplitudes arbi- 
trarily to Sj = T [ b 2 1 ~ J Sj i m ; b 2 b\~ J Sj_i >m ; ... ; 
where T is a known linear function of its arguments. The 
agreement between the matter amplitudes, Sj >m and the 
galaxy amplitudes Sj in the APM Survey indicates that 
the biasing parameters are b = b\ ~ 1 and bk — 0, for k > 1. 
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Otherwise, one would have to admit that there is a conspir- 
acy so that after the biasing transformation the resulting 
galaxy amplitudes Sj, at each order J, end up being similar 
to the matter ones. This is a strong constraint on bk even if 
the agreement were only within a factor of 2, e.g. the first 
term in T for J = 7 constraints 6 5 <2or6<1.15. 

To estimate Sj from the angular data, sj, we have as- 
sumed that absolute magnitudes and positions of galaxies 
are uncorrelated. Is this reasonable? In this paper we have 
also performed an identical analysis (to the one in the APM) 
using the angular positions of galaxies in the CfA and SSRS 
magnitude limited catalogs (§4). The results found for S3 
are very similar to the ones found for the same catalogs in 
redshift space (Gaztanaga, 1992). Moreover, they are also 
similar to the ones found in the angular distribution of vol- 
ume limited samples (Fry & Gaztanaga 1994) which, to first 
approximation, are not affected by redshift distortions or 
luminosity-density assumptions. The agreement between all 
these estimations indicates that luminosity-density correla- 
tions do not affect much the values of S3 or S4 obtained 
from magnitude limited catalogs. 

The values of S3 inferred from the angular data in the 
CfA/SSRS are smaller that the ones in the APM. As ar- 
gued in §2.5 it is unlikely that this discrepancy is caused by 
uncertainties in the luminosity function or its evolution, be- 
cause the inversion coefficients, rj, are quite insensitive to 
this. In the analysis over the Perseus-Pisces Redshift Sur- 
vey by Bonometto et al. (1993), the values of Q3 seem to 
decrease with increasing scales between 1 and 5 ft -1 Mpc. 
The corresponding variation for S3 is larger than the one in 
Figure 3, while at large scales S3 seems smaller than in the 
APM Survey and close to the values in the CfA/SSRS. We 
have shown that the small amplitudes of S3 in the CfA and 
SSRS catalogues correspond to unlikely fluctuations in the 
scatter of similar APM zones (Figure 6). Thus, our analysis 
indicates that, to estimate high-order correlation functions, 
the volume of a "fair sample" has to be much larger that 
the one in the combined CfA/SSRS catalogues or in the 
Perseus-Pisces Redshift Survey. The agreement of S3 in the 
APM (southern galactic cap) with S3 in the Lick (north- 
ern galactic cap) confirms this idea and indicates that the 
volume of the APM is probably close to a "fair sample". 

Following this sampling argument one might expect 
that the APM results should be closer to the results in the 
IRAS catalogs, which cover a bigger volume of space than 
the CfA or SSRS. But, in fact, Meiksin et al. (1992) and 
Bouchet et al. (1993) have found that S3 is even smaller 
in the IRAS distribution. The galaxy densities in cores of 
clusters determined from IRAS galaxies are systematically 
lower than those determined from optically selected galaxies 
(Strauss et al. 1992). This effect alone can reduce the values 
of S3 by a factor of two (Bouchet et al. 1993). Gaztanaga 
(1992) has argued that there are indications of a non-linear 
relation between optical and IRAS fluctuations. Thus, IRAS 
catalogs can not be directly compared with optical ones. 

That the skewness in the nearby optical catalogues 
are significantly lower than in the APM Survey, might in- 
dicate that the excess of big voids found (in a "sponge- 
like" topology) in our local galaxy distribution, i.e. within 
T> ~ 50 ft -1 Mpc, is not a representative feature of the Uni- 
verse. 



To summarise, we believe that the clustering properties 
found in the APM Galaxy Survey support the idea that 
large scale fluctuations in the universe grow according to 
the gravitational instability scenario and that galaxies trace 
matter at large scales. The values of the amplitudes found 
up J = 9 at small scales might be used to test non-linear 
growth and models for structure formation. 
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APPENDIX Al: CONNECTED MOMENTS 
AND DISCRETENESS 

To estimate the connected graphs, it is convenient to intro- 
duce the moment generating function M(t): 
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Figure Al. (a) Comparison of the two-point correlations with 
and without the discreteness correction. Squares and the contin- 
uous line correspond respectively to ~2 and j-i2 / N in the entire 
map. Circles and the dashed line correspond to ~2 (0) and ji2 1 N 
in a diluted sample, (b) Same as in (a) for the three-point corre- 
lations. 
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= E ^ =< e " >. 



so that: 
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Because of the discreteness /tj =< 8 J g > c N J are not 
good estimators oiUjN J unless N ^> 1. A better estimator 
of TojN is given in terms of kj using the Poisson model 
introduced in §2.1. In this model, the net effect of the dis- 
creteness is that the generating function of the discrete field, 
Mpoi SS on(t) is: 



The connected moments are then obtained from: 
d J 



MJ 



dt J 



log M(t) 



(A3) 



Up to J = 9 we find: 



M Polsson (t) = M{e* - 1) 



(A5) 



(see Peebles 1980, §33 or Gaztanaga & Yokoyama 1993) 
Using this relations we find: 
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(A6) 



where terms to the right of jij are the shot-noise correction. 

We have checked the accuracy of this shot-noise model 
by comparing the correlations in a given APM sample with 
the ones in a diluted version of the same sample, where 
only one out of twenty galaxies are selected at random. This 
is illustrated in Figure Al, where the connected moments 
< Sg > c = fij/N are quite different at small scales in the 
diluted sample and in the entire catalogue, whereas To j , ob- 
tained from kj in equations (A6), agree well for all scales. 
(The APM sample and the details of the analysis are de- 
scribed in section S3.) 



Cj as for a given J there are more sets of labels in the 
"snake" configuration. The resulting Cj are only slightly 
bigger than unity. For example for J = 3 we find numer- 
ically C 3 ~ 1.0031,1.0087 and 1.0166 for 7 = 1.4,1.7 and 
2.0. In our estimations of Cj we take a weighted average 
between the "star" and "snake" graphs. We have checked 
this procedure with the exact results up to J = 5. Because 
of the uncertainties in the observations a more accurate de- 
termination is not necessary. The values we find for 7 ~ 1.8 
are: C 3 ~ 1.01, C 4 ~ 1.03, C 5 ~ 1.05, C 6 ~ 1.07, C 7 ~ 1.10, 
C's ~ 1.13 and C 9 ~ 1.16. 

A very similar analysis follows for Bj, in three dimen- 
sions. For spherical cells, B3 has an analytical expression: 
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Gaztanaga & Yokoyama (1993). For J > 3 we take a 
weighted average between the "star" and "snake" graphs as 
before. For 7 ~ 1.8 we find B 3 ~ 1.03, 5 4 ~ 1.10, B 5 ~ 1.15, 
B e ~ 1.25, B 7 ~ 1.35, B s ~ 1.44 and B 9 ~ 1.58. 



APPENDIX A2: COEFICIENTS Cj AND Bj 

Even for the case of a power-law correlation {2 ~ , pa- 
rameters Bj in equation (21) have no simple analytical ex- 
pression because of the different topologies in the reassign- 
ments of labels a, b in equation (18). The same happens with 
Cj in equation (29) which are just the angular analogue of 
Bj with w 2 ~ e 1 "*. 

Two important classes of graphs can be estimated for 
any order. One is the "star" graphs where one of the labels 
in the pair is fixed, i.e. (1, 3)(2, 3)...( J, 3). The other is the 
"snake" graph where (1, 2)(2, 3)...( J - I, J). For the "star" 
graph we find: 



Cj(l) = -fe , F k = 2 zdzGtH*) (A7) 
^2 Jo 
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63(2) = I x dx I dcf> (z 2 + x 2 — 2zx cos (f>f 1 ~" t ^ 2 
Jo Jo 

which is all we need up to J = 3, i.e. C3 = C* and B3 = B*. 
For the "snake" graphs with J > 3: 

C'j(y) = F 2 ~ J+1 I zdzGj-^z) Gj(z) (A8) 
Jo 
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Gj(z) = I x dx I d(f> Gj-\(x) 
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/ 2 1 2 „ J.\(l-7)/2 

(z + x — 2zx cos <py " 

The contribution for J = 4 is C 4 = (12C| + 4C 4 )/16. For 
J > 4 more graphs have to be considered in a straight 
forward way. In general the numerical values for Cj are 



